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1.  Introduction 


This  report  considers  the  problem  of  distributed  detection  in  a  large,  random  network  of  sensor 
nodes.  We  consider  a  localized  event  and  focus  on  the  case  where  sensor  nodes  communicate 
binary  decisions  over  a  single-hop  wireless  network  to  a  common  fusion  node.  The  fusion  node 
combines  the  received  information  in  order  to  make  a  global  decision  on  the  presence  or  absence 
of  a  signal  source. 

This  work  is  related  to  and  draws  from  recent  works  in  distributed  detection.  Niu  et  al.  [1]  show 
that,  if  local  decisions  are  communicated  perfectly,  the  fusion  rule  for  independent  and  identically 
distributed  (i.i.d.)  binary  observations,  conditioned  on  the  true  hypothesis,  simplifies  to  counting 
the  number  of  detections.  In  reference  2,  Niu  et  al.  provide  analytic  and  approximate 
expressions  for  the  counting  rule  in  a  large,  random  sensor  network  and  a  random  target  location. 
Chang  et  al.  [3]  incorporate  a  random  access  protocol  for  a  fixed  number  of  sensor  nodes,  and 
they  derive  a  fusion  rule  that  is  a  weighted  sum  of  the  number  of  l’s  and  0’s  successfully  received 
at  the  fusion  node.  Similarly,  Kapnadak  et  al.  [4]  incorporate  a  random  access  protocol,  but  for  a 
random  sensor  network.  Aldalahmeh  et  al.  [5]  develop  a  real-time  counting  rule  and  provide  an 
approximation  to  the  system  performance  assuming  that  collisions  are  negligible.  Similar  to 
references  2,  4,  and  5,  we  consider  a  large,  random  sensor  network  and  incorporate  a  random 
access  protocol.  In  contrast,  however,  we  assume  a  sensor  node  transmits  only  when  it  declares  a 
detection  (i.e.,  when  a  node  decides  a  signal  source  is  present).  As  a  result,  collisions  from 
simultaneous  communications  imply  that  at  least  two  sensor  nodes  are  attempting  to  transmit 
detection  messages.  While  other  works  attempt  to  fuse  local  decisions  through  counts  of  l’s  and 
0’s  that  are  successfully  received,  this  work  derives  the  fusion  rule  for  successfully  transmitted  l’s 
and  a  count  statistic  on  collisions.  For  comparison,  we  re-derive  the  fusion  rule  assuming  the 
local  decisions  are  communicated  perfectly  [2]. 

The  remainder  of  this  report  is  outlined  as  follows.  In  section  2,  we  outline  the  sensor  network 
model,  including  the  local  sensor  operation,  and  the  network  communication  model.  Then  in 
section  3,  we  derive  the  fusion  rules  for  the  case  of  perfect  communications  and  a  given 
delay-constrained  communications  protocol.  In  section  4,  we  discuss  the  system  performance  as 
a  function  of  the  density  of  the  sensor  network.  We  derive  a  bound  on  the  vertical  uncertainty  of 
the  receiver  operating  characteristics  in  section  5.  We  present  numerical  examples  in  section  6 
validating  the  analytic  results  and  comparing  with  previous  results  under  perfect  communications. 
Finally,  conclusions  are  given  in  section  7. 
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2.  Network  System  Model 


In  this  section,  we  describe  the  sensor  network  model,  including  the  local  sensor  operation,  the 
network  communication  model,  and  the  fusion  node  processing.  We  assume  a  large  number  of 
sensor  nodes  are  randomly  deployed  independently  and  uniformly  over  a  bounded,  circular 
region.  Each  node  awakens  with  some  probability  at  each  time  interval,  and  if  active,  senses  its 
environment,  computes  a  local  decision,  and  communicates  any  detections  to  the  fusion  node 
using  a  delay-constrained  media  access  control  (MAC)  protocol.  When  multiple  nodes  attempt 
to  communicate  over  a  shared  wireless  communications  medium,  transmitted  messages  are 
subject  to  errors  due  to  collisions  with  other  ongoing  message  transmissions.  Therefore,  the 
fusion  node  must  take  this  into  account  when  forming  a  decision  rule. 

As  a  matter  of  notation,  unless  otherwise  specifically  stated,  random  variables  are  denoted  with 
uppercase  letters,  specific  outcomes  are  represented  by  lowercase  letters,  and  bold  letters 
represent  vectors.  For  example,  X  is  a  random  variable  while  X  =  [Xr., . . . ,  X„]T  is  a  random 
vector  with  exactly  n  elements. 

2.1  Sensor  Network  Model 

Sensor  nodes  are  deployed  randomly,  independently  and  uniformly  over  circular  region  TZ. 
Without  loss  of  generality,  we  assume  TZ  E  M2  is  centered  at  (0, 0)  with  given  radius  r  <  oo. 
Figure  1  depicts  a  random  sensor  network  deployment  with  a  fusion  node  centered  in  region  TZ. 
The  locations  of  the  signal  source  and  sensor  i  are  denoted  by  Ls  and  L  -,  respectively,  with 
Euclidean  distance  between  them  given  by  A  =  ||A  —  Lj|.  The  sensor  locations  {L^}  are  i.i.d. 
with  known  distribution  -  uniform  on  TZ.  It  follows  that  for  any  given  location  of  the  signal 
source,  the  distances  {A}  between  the  signal  source  and  the  sensors  are  conditionally  i.i.d. 

To  conserve  energy,  each  sensor  node  follows  a  random  sleep/wake  schedule,  independent  of  all 
other  nodes.  Fet  N0  =  {0, 1,2,.. .},  the  set  of  non-negative  integers.  The  number  of  active 
nodes  N  during  any  activity  period  is  modeled  according  to 

xn 

Pr  (N  =  n)  =  e~x— ,  n  E  N0,  (1) 

n\ 

where  A  E  (0,  oo)  is  the  average  number  of  active  nodes  in  region  TZ.  We  assume  activity  periods 
are  non- overlapping  and  have  a  fixed  and  known  duration.  From  the  above  properties,  active 
sensor  nodes  represent  a  homogeneous  Poisson  point  process  (PPP). 
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Figure  1 .  A  random  sensor  network  deployment  with  links  to 
a  fusion  node  centered  in  region  7 Z  with  radius  r. 


Each  active  sensor  node  senses  its  environment  and  makes  a  local  detection  decision.  Sensors 
collect  noisy  measurements  of  the  environment  for  the  purpose  of  detecting  whether  or  not  a 
source  signal  is  present  in  the  scene.  Under  the  null  hypothesis  ho,  when  no  signal  is  present, 
sensors  measure  only  noise.  We  assume  the  observations  are  i.i.d.  given  h0  is  the  true  state  of 
nature.  The  conditional  joint  distribution  of  the  sensor  observations  under  h0  is  given  by 

n 

fx\N,H  (x\n,  ho)  =  ]^[  fx\H  (Xi\h0) ,  (2) 

2—1 


where  X_  =  [A"i, . . . ,  Xn]T,  x  =  [xi, . . . ,  xn]T  with  ^61  for  7  =  1,2,...,  n,  and  x  is  a 
realization  of  X.  We  do  not  consider  the  more  general  case,  where  Xi  €  Rrn  with  m  >  1.  Here, 
Xi  could,  for  example,  represent  the  output  of  an  energy  detector  at  node  i. 

Under  the  alternative  hypothesis  h\,  in  which  a  source  is  present,  each  sensor  measures  a  source 
signal  embedded  in  noise.  We  model  the  observations  under  h\  to  be  conditionally  i.i.d. 
according  to 


n 

fx\A,N,H  ( x\a ,  71,  hi)  —  fx\A,H  (xi\a,  hi) ,  (3) 

2—1 


given  7i  active  nodes  and  the  signal  vector  a  =  [a  | , . . . ,  an]T,  where  a  is  a  realization  of  the 
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random  signal  vector  A  —  [Ai, . . , ,  /ln]T.  Each  A,  is  an  attenuated  version  of  the  scalar  signal 
level  As  >  0  originating  from  the  source.  We  assume  a  homogeneous  medium  that  attenuates  the 
source’s  signal  in  a  manner  that  does  not  depend  on  the  orientation  relative  to  the  source.  We 
model  the  signal  at  sensor  i  as  At  :=  Asg  (Z>0  >  0,  where  g  :  [0,  oo)  — >  (0,  oo)  models  the 
attenuation  as  the  signal  propagates  from  the  source  and  is  a  non-increasing  function  of  distance. 
The  attenuation  factor  g  can  model,  for  example,  spherical  spreading. 

We  assume  the  source’s  level  As  and  location  Ls  are  random  with  known  distributions  and  are 
mutually  independent  along  with  N,  { Lt } ,  and  sources  of  noise  in  sensor  observations.  Since 
{Lj}  are  i.i.d.,  it  follows  that  the  random  variables  Ai  are  identically  distributed,  while  { A, }  are 
not  necessarily  mutually  independent.  Let  ( Rs ,  ©s)  represent  the  polar  coordinates  of  the  signal 
source  located  at  Ls.  It  was  shown  in  reference  6  that  { At }  are  conditionally  independent  of  the 
angle  ©s  and  therefore  i.i.d.  given  Rs,  As,  and  N.  We  define  a  set  V  =  (0,  oo)  x  [0,  r]  and 
signal  parameter  vector  $  =  [/l,,  RS)T  e  V  with  probability  distribution  /$.  The  conditional 
joint  distribution  of  the  signal  at  the  sensors  is  then  given  by 

n 

fA\$,N,H  (a |0,  n,  hi)  =  ]^[  }a\®,h  (aj|0,  hi)  .  (4) 

i= 1 

The  vector  $  captures  the  signal  source  components  that  are  shared  through  observations  at  all 
active  sensor  nodes.  One  consequence  is  that  such  a  network  could  infer  only  the  radial 
component  of  the  signal  source’s  location  from  sensor  observations.  In  a  practical  setting,  the 
sensor  deployments  and  signal  propagation  loss  functions  may  not  be  rotation-independent.  In 
which  case,  orientation  may  also  be  inferred  from  observations  given  an  appropriate  model. 

While  $  (and  realization  0)  represents  a  vector  of  the  signal  source  parameters,  specifically  the 
signal  level  and  range  to  the  origin  of  7 Z,  we  loosely  refer  to  $  (and  0)  as  the  signal. 

For  the  given  model,  it  follows  that  the  sensor  observations  are  conditionally  i.i.d.  according  to 

n 

Sx\*,n,h  (z|0,  n,  hi)  =  JJ  fx\*,H  (a?t|0,  h)  ,  (5) 

i=l 

where  fX\*,H  (®|0,  ^i)  =  E  (fX\A,  (x\A,  0,  h\))  and  E(-)  denotes  expectation  (it  is  with 
respect  to  zl|$,  H  in  this  case).  The  probability  distribution  fx\$,H  (a’;|0,  hi)  corresponds  to  the 
conditional  distribution  of  an  observation  at  active  node  i  under  hi  and  given  the  signal  0.  This 
sensor-level  distribution  is  essentially  the  result  of  marginalizing  the  randomness  in  a  sensor 
node’s  location  -  the  randomness  that  remains  in  Ai  after  conditioning  by  $.  Note  that  { X, }  are 
conditionally  independent  of  $  given  h0  is  the  true  state.  Thus,  fx\&,H  (x\4>,  h0)  =  fx\H  (x\h0). 
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It  follows  that  under  both  hypotheses,  the  observations  are  conditionally  i.i.d.  across  the  sensor 
network  for  a  given  signal  0  and  n  sensor  nodes. 

Each  active  sensor  node  first  makes  a  local  binary  decision  and  then  communicates  that  decision 
to  a  common  fusion  node.  A  local  decision  is  a  quantized  or  compressed  version  of  the  of  the 
original  observation.  We  consider  two  communications  modes:  a  perfect  communications 
model  and  a  delay-constrained  random  access  communications  protocol. 

2.2  Distributed  Detections 

Each  active  sensor  node  makes  a  local  binary  decision  and  communications  that  decision  to  a 
common  fusion  node.  The  local  decisions  are  determined  by  an  ctj-level  decision  rule,  denoted 
by  The  local  decision  rules  d,  are  not  necessarily  identical  across  the  sensor  network  in 
optimal  fusion  of  binary  decisions,  even  when  observations  are  conditionally  i.i.d.  [7].  However, 
we  assume  each  sensor  shares  a  common  o-levcl  rule  S.  This  choice  dramatically  simplifies  the 
design  and  deployment  of  the  sensor  network,  and  in  some  cases,  proves  to  be  asymptotically 
optimal  for  binary  decisions  [8]. 

The  parameter  a  £  (0,1)  represents  the  local  probability  of  false  alarm.  The  local  probability  of 
missed  detection,  denoted  by  /3,  is  given  by 

/3  =  l~E(6(X)\h1).  (6) 

We  denote  the  local  decisions  by  Yt  £  {h(j.  h{\  for  i  =  1,  2, . . . ,  N.  We  refer  to  Y]  =  hi  as  a 
local  detection,  whether  it  is  a  true  (H  =  hi)  detection  or  a  false  (H  =  h0)  detection.  Each  Yt 
conditioned  on  <3>,  N,  H  can  be  seen  as  a  Bernoulli  random  variable  with  success  probability 
Pj  (S  |0)  defined  by 


pj(6\4)=E(6(X)\l,hj),  (7) 

for  j  £  {0, 1}.  The  success  probabilities  can  be  viewed  as  the  conditional  probabilities  of  false 
alarm  (j  =  0)  and  detection  (j  =  1)  for  decision  rule  5  and  signal  cf).  When  convenient,  we 
suppress  the  dependence  on  5  by  the  local  probabilities  of  detection  and  false  alarm  and  simply 
write  pi  (0)  =  pi  (5 10)  and  p0  (0)  =  p0  (<5|0) ,  respectively.  Note  that  value  p0  (0)  =  a  and 
does  not  depend  on  0,  while  the  decision  region  for  S  may. 
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The  conditional  joint  distribution  of  the  binary  decisions  Y_ |$,  N,  H  is  then  given  by 

W,HW,n,^  =  (^r[1_ftWr"’  **"  (8) 

10,  z  >  n, 

for  n  £  N0,  j  £  {0, 1},  and  where  ;j  =  Y^=i  1  (Vi  =  hi)  and  l(-)  is  the  indicator  function.  The 
joint  probability  distribution  of  Y,  7Vj0,  H  is  then 

_  \n 

fv,N\YH (y,  n\<j>,  hj )  =  Pj  (0) *  [1  -  Pj  (0)]  n~z  e~A— ,  (9) 

for  z  <  n.  Let  Z  =  Y)f=i  1  {Yi  =  hi).  For  n  active  nodes  and  z  <  n  detections,  there  are  (") 
different  y  that  satisfy  J0"=i  l(yt  =  h  \  )  =  z  with  equal  probability.  It  follows  that  the  probability 
distribution  of  Z,  N\H  is  given  by 

fz^H{z,n\^hj)  =  (^JPj  &  tf)r~*e~A“j  »  (10) 


for  z  <  n  <  oo.  Marginalizing  N,  the  number  of  active  nodes,  the  conditional  distribution  of 
Z|$,  H  is  then  given  by 


/z|$„ff(0|0,  hj) 


_  c-Ap,-fo)  \-Xpi  Ml 

z\ 


(11) 


for  z  £  No-  Note  that  the  dependence  on  $  in  equations  7  through  1 1  drops  out  under  hypothesis 

h0. 


Recall  that  the  active  sensor  network  represents  a  homogeneous  PPP  and  the  number  of  active 
nodes  is  Poisson  with  mean  A.  It  follows  by  the  thinning  property  of  a  PPP  that  Z |$,  H  is  a 
Poisson  random  variable  with  parameter  A  pj  (0)  under  hj  and  given  $  =  0. 

2.3  Network  Communication:  Protocol  Model 


In  this  section,  we  describe  a  delay-constrained  MAC  protocol  model.  The  local  decisions  Yt,  for 
i  —  1,  2, , . . ,  N,  are  communicated  to  a  fusion  node.  The  fusion  node  then  combines  received 
information  to  make  a  global  decision  as  to  the  presence  or  absence  of  a  signal  source.  The 
global  decision  rule  must  account  for  the  effects  in  communications. 

Sensor  nodes  are  either  awake  or  asleep  during  an  activity  period,  which  is  formed  by  a  sensing 
period  followed  by  a  communicating  period.  The  sensing  and  communicating  periods  are 
synchronized  so  that  each  active  sensor  node  observes  the  environment  and  makes  a  local 
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Figure  2.  Activity  periods  composed  of  a  sensing  period  and  a  communicating 
period  for  all  active  sensor  nodes. 


decision  during  the  sensing  period.  Then,  nodes  attempt  to  transmit  their  decisions  to  a  common 
fusion  node  during  the  communicating  period.  Figure  2  graphically  depicts  the  sensing  and 
communicating  periods.  Note  that  there  could  be  overlap  between  a  communicating  period  and 
the  next  sensing  period,  assuming  one  does  not  interfere  with  the  other  and  so  that  similar  periods 
do  not  overlap  (e.g.,  Sensing  Period  t  does  not  overlap  with  Sensing  Period  t  +  1).  We  model  the 
system  as  being  time-invariant.  Thus,  we  need  only  model  a  single  activity  period. 

To  both  save  transmission  energy  and  reduce  communication  collisions,  sensor  nodes  only 
transmit  detections  (i.e.,  when  Yt  =  hi).  Furthermore,  the  detection  messages  do  not  need  to  be 
unique  to  the  sensor  nodes.  In  references  4  and  5,  similar  distributed  detection  models  were 
considered,  but  a  transmitted  decision  can  represent  either  h0  or  h , .  Here,  messages  to  be 
transmitted  represent  only  Yt  =  h\  decisions.  In  addition,  we  assume  sensor  nodes  do  not 
retransmit  messages  in  the  event  of  a  collision,  again  to  conserve  energy.  Since  transmitted 
messages  only  represent  detections,  collisions  provide  valuable  information  at  the  fusion  node. 

During  the  communicating  period,  each  sensor  with  a  detection  attempts  to  transmit  a  single 
detection  message  via  Slotted  ALOHA  (S-ALOHA).  The  communicating  period  is  broken  into 
1  <  M  <  oo  equal-duration  slots.  On  the  receiving  side,  the  fusion  node  does  not  know  which 
slots  will  be  used,  so  it  must  wait  the  duration  of  all  M  slots  during  the  communicating  period 
before  making  a  decision.  Thus,  the  communications  delay  is  directly  proportional  to  M .  The 
delay,  and  therefore  the  number  of  slots  M,  is  specified  according  to  the  needs  of  the  detection 
application.  In  a  slight  deviation  from  previous  notation,  M  is  not  considered  random. 

Nodes  with  detections  randomly  choose  one  out  of  M  slots,  independent  of  all  other  nodes, 
according  to  the  same  probability  distribution.  Let  Si  represent  the  slot  number  chosen  by  the  ith 
node.  The  slot  numbers  are  then  drawn  according  to  Pr  (S'*  =  m)  =  pm,  for  rri  —  1,2, ... ,  M. 

We  assume  the  fusion  node  can  detect  collisions  in  any  of  the  M  slots.  A  message  is  considered 
to  be  received  successfully  in  a  slot  if  only  one  sensor  node  transmits  in  that  slot.  Otherwise,  a 
slot  is  either  unused  or  contains  collisions.  A  slot  with  2  or  more  transmissions  is  counted  as  1 
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collision.  We  denote  the  number  of  sensor  nodes  that  select  slot  m  for  transmitting  a  detection 
message  by  Tm,  and  Um  represents  the  observation  in  slot  m  at  the  fusion  node.  We  also  refer  to 
Tm  as  the  slot  occupancy  (i.e.,  the  number  of  nodes  attempting  a  transmission  in  the  slot).  The 
observations  at  the  fusion  node  are  then  related  to  the  slot  occupancies  according  to: 


T  0,  if  Tm  =  0 

Um=l  1,  if  Tm  —  1  (12) 

[  c,  if  Tm  >  1, 

for  m  —  1,  2, . . . ,  M  and  where  c  signifies  a  collision  in  the  slot.  With  regards  to  notation,  the 
vector  U_  =  [U\,  f/2, . . . ,  Um]t  represents  the  collection  of  observations  at  the  fusion  node  at  the 
end  of  the  communicating  period,  and  u  =  [ui,  u2, . . . ,  um]t  is  a  particular  realization  of 

U  G  (0, 1,  c}M. 


Given  the  number  of  nodes  with  transmissions,  the  slot  occupancies  do  not  depend  on  the  signal 
$  or  the  hypothesis  H.  Since  the  slot  numbers  are  selected  i.i.d.  across  the  nodes  with 
detections,  the  probability  of  occupancy  in  slot  m,  conditioned  on  z  transmissions,  is  given  by 


hm\z  (k\z) 


(t)i( 1  -  PmY  k,  k  <  Z 
0,  k  >  z, 


(13) 


for  k,z  G  No.  Under  this  network  model,  the  system  designer  can  select  the  distribution  from 
which  nodes  choose  slot  numbers.  For  a  discrete  uniform  distribution  (i.e.,  pm  =  1/M  for 
m  —  1,2,...,  M ),  the  slot  occupancies  become  identically  distributed  across  slots  and  the 
number  of  transmitting  nodes  are  uniformly  spread  across  slots.  We  assume  the  slots  are  chosen 
according  to  the  uniform  distribution  for  the  remainder  of  this  work. 


The  number  of  active  nodes  is  unknown  to  the  fusion  node  and  the  number  of  transmitting  nodes 
are  uniformly  spread  across  slots.  The  slot  occupancies  represent  a  thinning  of  Poisson  points 
represented  by  sensor  nodes  with  detections  by  the  proportion  jj-  It  is  straightforward  to  show 
that  the  slot  occupancies  are  conditionally  i.i.d.  and  distributed  according  to 


fTm\^H  {k |0,  hj) 


-  c-TTP:,M  [MPj  Ml 
k\ 


(14) 


for  k  G  N0.  The  fusion  node  observes  U_  G  (0, 1,  c}M  with  slot  occupancies  mapped  according  to 
equation  12.  The  probability  distribution  of  U_ |$,  H,  after  some  simplification,  is  then  given  by 


fu\$,H  (m|0,  hj) 


(^Oj 


\M—n\—nc 


(TTij)”1  (7T, 


c,JJ 


(15) 
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where  n \  =  Yhm=i  l(wm  =  1)  and  nc  =  Y^m=i  1 (um  —  c )  are  just  the  counts  of  I  s  and  c’s  across 
slots,  respectively,  and 

7T0j  =  Pr  (Tm  =  0|  0,  hj ) 

7Tij  =  Pr  (Tm  = 

7rCii  =  Pr  (Tm  >  1|0,  hj)  .  (16) 

The  probability  nuj  depends  on  the  local  decision  rule  5  under  both  hypotheses  and  the  signal  $ 
under  h  i .  These  dependencies  are  not  explicit  in  denoting  the  probabilities  nuj  to  simplify  the 
notation,  but  are  clear  from  equation  14. 

We  define  the  (random)  counts  Nu  =  1(^4 n  —  u).  It  is  clear  from  equation  15  that  the 

observation  {N\,  Nc}  is  equivalent  to  U_  for  determining  H .  The  observations  at  the  fusion  node 
are  then  just  the  counts  of  the  I  s  and  c  s. 


3.  Global  Decision  Rules 


In  this  section,  we  derive  the  global  decision  rule  (or  fusion  rule)  for  a  given  local  decision  rule. 
First,  we  derive  the  fusion  rule  assuming  the  local  decisions  are  perfectly  communicated  to  the 
fusion  node.  Then,  we  derive  the  fusion  rule  for  the  delay-constrained  MAC  protocol  outlined  in 
section  2.3. 

The  global  decisions  rules  are  determined  by  an  o-lcvcl  composite  likelihood-ratio  test.  The 
form  of  the  rule  is  the  same  throughout,  given  by 

(  1,  if  l(i)  ( B )  >  7 7(i) 

%)(£)=]  C(i),  if  l(i)  (5)  =  77(i)  (17) 

{  0,  if  l(i)  (B)  < 

where  ( i )  denotes  the  special  case,  B  is  drawn  according  to  a  probability  distribution  fs\H,  and 
/(q  is  the  likelihood  ratio.  The  probabilities  of  false  alarm  and  missed  detection,  respectively,  are 
given  by 


«(*)  =  E  (%)(£)  | /i0) 


(18) 
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and 


/%)  =  1  —  E  (S^(B)\hi)  .  (19) 

The  threshold  rj^  >  0  and  randomization  parameter  Qt)  e  [0, 1]  are  chosen  to  satisfy 

min/5(j)  s.t.  a (q  =  a,  (20) 

5(i) 

for  some  a  £  (0, 1).  All  fusion  rules  are  subject  to  the  same  global  false  alarm  probability  a. 

The  global  rules  are  potentially  suboptimal  since  the  local  decision  rules  are  constrained  to  be 
identical  across  the  network. 

3.1  Perfect  Channel 

We  first  consider  fusing  the  binary  decisions  from  all  active  sensors  assuming  unconstrained  and 
error-free  communications.  In  other  words,  it  is  assumed  the  local  decisions  are  received 
perfectly  at  the  fusion  node. 

The  fusion  node  observes  Z  £  N0,  the  number  of  active  nodes  with  detections.  From  equation 
11,  the  conditional  likelihood  ratio  for  the  global  decision  rule  5(i)  is  given  by 

l{1)(z\D  =e-AMtM 

for  a,a  £  (0, 1).  The  threshold  rj (i)  and  the  randomization  parameter  Q , )  are  chosen  to  satisfy 
the  d-level  test  with  likelihood  ratio  l(i)(z)  =  E(((i)(z|$)). 

3.2  Delay- Constrained  MAC 

Now  we  consider  the  fusion  rule  for  the  delay-constrained  MAC  protocol.  The  binary  decisions 
are  communicated  to  the  fusion  node  with  sensor  nodes  transmitting  according  to  the  MAC 
protocol  over  M  slots.  The  fusion  node  forms  the  decision  rule  after  receiving  the  observations 
upon  completion  of  all  M  slots. 

The  fusion  node  observes  iVi,  iVc  e  (0, 1, . . . ,  M}  with  Ni  +  Nc<  M.  The  fusion  node  also 
observers  the  number  of  unused  slots  N0,  which  is  a  known  quantity  given  Nu  Nc,  and  M.  From 


Pi 


a 


(21) 
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equation  15,  the  conditional  likelihood  ratio  for  the  global  decision  rule  <5(2)  is  given  by 


k 2) 


7T0,1 

M—n\—nc 

’tt  i.r 

ni 

7Tc,l 

.7^0,0. 

_7Tl,0_ 

_^"c,  0 . 

(22) 


for  ri\  +  nc  <  M,  ttuj  /0,a6  (0, 1)  and  a  G  (0, 1).  The  threshold  r/( 2)  and  the  randomization 
parameter  Q2)  are  chosen  to  satisfy  the  a-lcvcl  test  with  likelihood  ratio 

1(2  )(ni,nc)  =  E(/(2)(ni,nc|$)). 


The  likelihood  ratios  for  both  cases  are  summarized  in  table  1 .  The  expectations  are  with  respect 
to  the  random  signal  $.  In  each  case,  the  decision  rules  are  functions  of  count  statistics:  the 
number  of  active  nodes  with  detections,  the  number  of  slots  with  a  successful  transmission,  or  the 
number  of  slots  with  collisions. 


Table  1.  Likelihood  ratios  of  count  statistics. 


Perfect  Communications 

MAC  Protocol 

ki)(z)=E\ 

^e-A[pi($)-a] 

Pi® 

a 

') 

i 

l(2)(ni ,nc)  =  E  | 

r 

7T0,1 

_7r°,°_ 

M—n\—nc 

7T1,1 

7Tl,0 

ni 
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4.  System  Performance 


In  this  section,  we  discuss  the  system  performance  as  a  function  of  the  density  of  the  sensor 
network  and  the  number  of  communications  slots.  The  coverage  area  remains  fixed,  thus  the 
average  number  of  active  sensors  is  synonymous  with  density.  The  system  performance  is 
characterized  in  terms  of  average  energy  consumption  and  global  probability  of  missed  detections. 

4.1  Energy  Versus  Sensor  Network  Density 

The  amount  of  energy  consumed  for  wireless  transmissions  can  be  of  the  same  order  as  that 
needed  for  sensing  and  processing  [9].  For  the  analysis  here,  we  simply  assume  these  two  levels 
are  identical.  A  more  thorough  analytical  study  or  one  like  reference  9  could  be  carried  out  for 
the  specific  system  being  deployed. 

Let  E  represent  the  amount  of  energy  consumed  by  a  sensor  node  over  a  single  sensing  or 
communicating  period.  Let  p0  =  Pr  (H  =  h0).  It  is  straightforward  to  show  that  the  average 
energy  consumption  across  the  sensor  nodes,  denoted  Es,  is 

Es  =  A[1  +  poa  +  (1  —  p0)(l  —  P)\E,  (23) 

where  a  and  (3  are  the  local  probabilities  of  false  alarm  and  missed  detection,  respectively. 

We  assume  that  E  is  constant  with  sensor  network  density.  In  the  single-hop  setting  here,  the 
coverage  area  and  the  location  of  fusion  node  are  fixed.  It  is  assumed  that  E  is  the  lowest  amount 
needed  so  that  any  one  node  in  the  coverage  area  can  reliably  communicate  with  the  fusion  node. 
Collisions  are  permitted  and  no  attempt  is  made  to  separate  near  and  far  transmitters.  A  constant 
E  could  have  an  impact  on  the  fusion  node’s  ability  to  differentiate  between  the  three  states  in  a 
slot:  no  transmission,  one  transmission,  or  a  collision.  We  do  not  consider  physical  layer 
complexities  here. 

The  fusion  node  must  remain  active  for  all  communicating  periods.  The  number  of  slots  and 
communications  rates  remain  fixed  with  respect  to  density.  Thus,  the  amount  of  energy 
consumed  by  the  fusion  node  remains  constant  with  density.  Let  Ef  denote  the  amount  of  energy 
consumed  by  the  fusion  node  in  any  given  activity  period.  The  total  average  energy  consumed  by 
the  sensor  network  in  any  particular  activity  period  is  then  given  by 

Et  —  Es  +  Ef.  (24) 
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We  assume  the  signal  source  is  mostly  absent  (i.e.,  7 r0  >  1/2).  Most  activity  periods  do  not 
contain  a  signal  source,  and  so  processing  and  message  transmissions  are  usually  due  to  false 
alarms.  We  constrain  the  sensor  network  by  forcing  Aa,  the  average  number  of  false  alarms,  to 
be  constant  as  the  density  of  the  sensor  network  grows  (i.e.,  a  — >  0  as  A  — >  00  s.t.  Aa  is 
constant).  In  a  typical  setting,  /3  tends  to  1  as  a  tends  to  0.  It  follows  then  that  the  average 
per-node  energy  use  declines  as  the  sensor  network  density  increases.  So  while  the  total  average 
energy  across  the  network  grows  with  density,  the  per-sensor-node  usage  decreases  and  the  fusion 
node’s  usage  is  fixed. 

4.2  Error  Probability  Versus  Sensor  Network  Density 

Here,  we  analyze  the  global  probability  of  missed  detection  while  increasing  sensor  network 
density  and  fixing  the  average  number  of  false  alarms  from  the  sensor  nodes  and  fixing  the  global 
probability  of  false  alarm.  The  probability  of  missed  detection  is  discussed  for  both  cases. 

Condition  1.  The  local  detection  tests  are  such  that  the  local  conditional  probability  of  detection 
Pi  (5|0)  is  concave,  non-decreasing  in  a  G  (0, 1)  for  all  signals  0  G  V. 

Condition  2.  The  sensor  observations  under  h0  and  h\  have  the  same  support  for  all  0  V.  This 
condition  ensures  that  the  limit  limQ^0  px  (0)  exists.  The  following  results  can  be  easily 
extended,  however  tedious,  with  this  condition  removed. 

4.2.1  Perfect  Channel  with  Observation  { Z } 

Since  p\  (0)  >  a,  V0  G  V,  the  conditional  likelihood  ratio  l^)(Z\(f>)  is  non-decreasing  in  Z.  An 
equivalent  test  for  this  case  is  given  by 


11,  if  Z  >  z(  1) 

h)(Z  10)  =  <  C(i)i  ifZ  =  z(1)  (25) 

[0,  ifZ<  z( i), 

where  now  zm  and  are  chosen  to  satisfy  d.  It  is  clear  that  the  decision  region  for  this  d-level 
test  does  not  depend  on  0.  Thus,  the  decision  region  for  rule  <5(i)  (Z|0)  that  achieves  d  under  ho 
does  so  also  for  composite  rule  0i) (Z).  Since  the  decision  region  for  5(!) (Z|0)  does  not  depend 
on  0,  we  have  (3(  1)  =  E (0(!)($)). 
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Let  0  €  V  be  given.  The  conditional  probability  of  missed  detection  is  given  by 


2(o 

Ai)«  =  Ee~'Vl(£l 


[Api  (0)  ' 


2  =  0 


Z! 


-  C(i)e 


-Api  (^) 


[A?i  (0)  ] 

-(i)! 


(26) 


If  o;  is  fixed  and  A  a  is  held  constant  as  A  increases,  then  from  equations  1 1  and  25  both  Z(i}  and 
C(1)  are  fixed.  Then,  any  change  in  ftp)  (0)  is  due  to  a  change  in  Xpi  (0) . 

Let  p  =  Xp\  (0) .  It  is  straightforward  to  show  that  equation  26  is  decreasing  in  //.  Thus,  if  p  is 
non-decreasing  in  A  for  every  0  6?,  then  /0|,  is  non-increasing  for  A.  It  remains  to  show  that 
the  average  number  of  detections,  Xp\  (0) ,  is  increasing  in  A  with  the  constraint  that  A  a:  remains 
fixed. 


Let  0  <  Ai  <  A2  <  00  and  0  <  a2  <  ot-i  <  1  be  such  that  Aiay  =  A2a2.  Additionally,  let  5i  and 
02  be  the  local  detection  rules  associated  with  a  1  and  a2,  respectively.  Since  the  local  tests  are 
concave  for  any  given  0  6?,  we  have 

?l(^l|0)  ^  ?1  (02 10)  /nn 

—  1  \A  • 

Q^l  Oi2 

with  equality  if  and  only  if  pi(<5|0)  is  linear  on  (0,  «i).  Subsequently,  we  have 


A2?i(02|0)  -  AiPi((5i|0)  >  [A2a2  -  Ai«i] 


Pi (^1 |0) 


=  0. 


(28) 


In  general,  the  average  number  of  positive  detections  is  non-decreasing  in  A  while  Ao  is  constant. 
The  average  number  of  positive  detections  is  constant  only  where  the  local  ROC  curves  have  a 
line  segment  meeting  the  origin.  If,  for  example,  the  local  ROC  curve  is  differentiable,  then  the 
average  number  of  positive  detections  is  strictly  increasing,  and  so  the  global  probability  of 
missed  detection  is  strictly  decreasing  in  A  for  a  fixed  average  number  of  false  alarms. 

4.2.2  MAC  Protocol  with  Observation  (iVi,  A0} 

Now,  we  analyze  the  performance  of  the  MAC  protocol  while  increasing  sensor  network  density. 

Let  p  =  Xpi  (0)  and  M  be  finite  and  fixed.  We  show  that  0(2)  is  decreasing  in  p.  The  conditional 
likelihood  ratio  /(2)(i,  j\4>)  is  strictly  increasing  in  each  coordinate  and  thus  l(2)('i-j)  is  strictly 
increasing  in  each  coordinate.  Also,  the  joint  probability  of  {A, ,  Nc  \  h() }  remains  fixed  with 
respect  to  p  and  does  not  depend  on  0.  Hence,  the  o-lcvcl  decision  region  does  not  depend  on  p 
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and  (j).  Hence,  it  suffices  to  show  that  /3(2)(0)  is  decreasing  in  /i. 

Let  X  =  {(i,j)  G  {0, 1, ,  M}2|i  +  j  <  M}  denote  the  feasible  set  of  points.  Also,  let 
A  =  { {%,  j)  G  X\l(2){i.j)  <  77(2) }  denote  the  a-lcvcl  acceptance  region  for  h0.  Under  the  MAC 
protocol,  the  global  probability  of  false  alarm  is  given  by 

M  M 

«(2)  =  J]  Ylb)  (Li)Pr(^!  =  i,Nc  =  j\  h0).  (29) 

i=0  )=0 

The  threshold  rjm  and  randomization  Q2)  satisfy  the  constraint  df(2)  =  a.  The  global  probability 
of  missed  detection,  conditioned  on  cf),  is  given  by 


M  M 


P{2)(<j>)  =  EE  (X  _  5(2)(L2'))  Pr(^i  =  i,Nc  =  j\hi,4>), 


i= 0  7=0 


which  is  upper  bounded  by 


M  M 

hM  <  E  EPr(7Vi  =  Nc = j \hi,t) 

i=o  j=o 
(iJ)eA 
M  M 

1 1  //  r 

(  \i(_  \j 
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i=0  7=0  v  J  J 
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M  M 
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^ ^ IV I  : 
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Denote  the  upper  bound  by  7i.  Differentiating  /i  with  respect  to  //,  we  have 


(30) 


(31) 


(32) 


(33) 


M  M 


^  EEPrW  =^,  =  7^,0)  (-1  +  1+  j  eM/M“1 

«=o  7=o  ' 
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d/x 


fi  Me»/M  -1-  n/M 


(34) 


The  joint  probability  masses  are  strictly  positive  since  7rU)i  >  0,  for  u  G  (0, 1,  c}.  Assume 
(0,  M)  cj  A  (Note:  This  point  maximizes  the  likelihood  ratio.  For  sufficiently  small  a,  we  can 
exclude  this  point  from  the  acceptance  region  for  h0 ).  Then,  for  a  weak  constraint  on  /i,  we  want 
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each  term  in  the  summation  to  be  negative.  Accordingly,  we  have  the  inequality 


n  1  |  i  |  j  —  1 

0>  ~1  +  Jl  +  Me#  -  1  -  p/M' 


(35) 


It  can  be  shown  that  the  second  term  on  the  right-hand  side  in  the  above  inequality  is  less  than  the 
third  term.  Thus,  any  p  that  satisfies  the  inequality  with  i  =  1  and  j  =  M  —  1  also  satisfies  the 
inequality  for  any  (i,j)  €  A.  After  rearranging  terms  and  weakening  the  above  inequality,  we 
have  the  inequality 


e^/M  -  2p  +  1  >  0,  (36) 

which  is  satisfied  for  sufficiently  large  p  since  lim,..^  e~xxK  =  0  for  any  fixed  and  finite  K .  In 
other  words,  the  global  probability  of  missed  detection  is  strictly  decreasing  in  p  provided  the 
probability  of  a  successful  transmission  in  any  slot  satisfies  7Ti.i  =  e~^!M (/x/ M)  <  Aj.  Note 
that  7Ti  i  is  strictly  decreasing  in  p  for  p  >  M. 

It  remains  to  show  that  the  point  (0,  M)  maximizes  the  likelihood  ratio.  It  is  sufficient  to  work 
with  the  conditional  likelihood  ratio  for  the  same  reasons  discussed  above.  Since  the  conditional 
likelihood  ratio  is  strictly  increasing  in  its  coordinates,  one  could  compare  each  of  the  M  +  1 
points  that  satisfy  ri\  +  nc  =  M  (i.e.,  the  upper  right  boundary  of  feasible  points).  Instead,  from 
the  likelihood  ratio  in  equation  22,  we  only  need  to  show  that  the  ratio  is 

greater  than  unity.  Assume  pi{4>)  >  a  for  at  least  one  (p  £  V  with  non-zero  probability. 
Otherwise,  pi((j>)  =  a  for  all  0  e  V,  with  all  the  local  tests  being  purely  random  guesses  (i.e.,  no 
sensor  data  required). 

Consider  that  0.  The  ratio  from  the  terms  of  the  likelihood  ratio  can  be  written  as 

ttc, 1*1,0  Zr-2 

TTc.oTTM  JX2  ^e-^/M(APl/M) 

pXpi/M  (Xpi/M)m 

Z—/m=l  (m+1)! 

pXa/M  V°°  (-WM)m 
£—/m=l  (m+1)! 

>  1, 

where  the  inequality  follows  since  f(x)  =  ex  i  (k+iy.  1S  clearly  increasing  in  x  e  R+.  Thus, 
the  point  (0,  M)  is  the  unique  feasible  point  that  maximizes  l(2)(ni,nc)- 


(37) 

(38) 

(39) 
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4.3  Error  Probability  Versus  Number  of  Communications  Slots 


In  this  section,  we  analyze  the  performance  of  the  fusion  rule  under  the  MAC  protocol  while 
increasing  the  number  of  communications  slots  M .  Instead  of  using  the  analytic  form  of  the 
probability  of  missed  detection,  we  examine  the  asymptotic  error  rate  for  a  given  signal  (j)  and 
fixed  and  finite  A.  For  that,  we  appeal  to  the  Chernoff-Stein  Lemma  [10]. 


For  a  given  <fi  e  V,  each  Ui,  U2, . . . ,  Um  are  drawn  i.i.d.  from  Q 0  or  Q i,  where  Qj  represents  the 
trinomial  distribution  of  U,  e  {0, 1,  c}  under  hj  for  j  =  0,1  and  z  =  0.1,...,  M.  According  to 
the  Chernoff-Stein  Lemma,  the  best  exponent  of  probability  of  missed  detection  is  given  by 

D(Qo\\Qi)  (40) 

^2  TTn.olog—  •  (41) 


^Mhe^=- 


Let  x0  =  Xa/M  and  X\  =  Xpi/M.  From  the  log-sum  inequality,  we  have  the  following  relations. 
The  relative  entropy  of  the  trinomial  is  lower  bounded  by 


o  XQ 


(1 


?—%o' 


D(Qo\\Qi)  >  e~x°  log  —  +  (1  -  e-*°)  log  (1  _  e_xi) 


OO  f, 

=  e~x°  |  xx  -  x0  +  (  x0  +  ^  I  log 

k= 2 


o  XO 


+ Er=2 1 


^  Z^fc= 2  k\ 


The  relative  entropy  of  the  trinomial  is  also  upper  bound  according  to 

OO  fc 

ZWolIQi)  <  5>-“fpog 


k= 0 


p  -Xo  ±0 

e  k\ 


p-x  1  ±X 
C  k\ 
X0 


=  Xi  ~  X0  +  X0  log  — 

X\ 

=  bDWW' 


(42) 

(43) 


(44) 

(45) 

(46) 


where  Pq  and  l\  are  Poisson  distributions  with  mean  A  a  and  Xp\  respectively.  From  equations 
43  and  46,  we  have 


lim  MD{Qq\\Q1)  =  D{Pq\\P1).  (47) 

M->  oo 

This  limit  is  not  unlike  the  performance  of  the  perfect  communications  case.  To  make  the  link, 
we  compare  the  asymptotic  behavior  of  the  MAC  protocol  to  that  of  another  system  under  perfect 
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communications  with  observations  that  converge  in  distribution  to  the  Poisson  case.  Namely,  the 
observations  seen  by  the  fusion  center  are  binomial. 


Let  n  be  the  number  of  sensor  nodes,  and  the  fusion  node  observes  a  sequence  Y\,  Y2, . . . ,  Yn. 
Each  Yi  are  conditionally  i.i.d.  and  drawn  from  either  B0  under  h0  or  B\  under  h\.  The 
distributions  B0  and  B\  correspond  to  binomial  distributions  with  parameters  ( n,paa )  under  ho 
and  (n,  papi),  where  pa  G  (0, 1)  is  the  probability  of  a  sensor  node  being  active.  We  let  pa  — >•  0 
as  n  — »  og  such  that  npa  =  A.  It  can  be  shown  that  the  statistic  Z  =  JA  Yt  is  a  sufficient  statistic 
for  detection  at  the  fusion  node.  It  is  well  known  that  the  conditional  distributions  of  Z  converge 
to  Poisson  distributions  with  means  A  a  and  Xp\,  respectively.  Thus,  the  performance  of  this 
system  converges  to  that  of  the  perfect  channel  system  discussed  in  previous  sections. 
Additionally,  the  best  exponent  of  the  probability  of  missed  detection  as  n  — >  oo  is  the  negative  of 
relative  entropy  given  by 


D(B0\\B1)  =  paa  log  +  (1  -paa)\og  1 


PaOi 


PaP  1 

=  paa  log  —  +  (1  -  paoi)  (  y ^ 


1  -  PaPl 
(PaPl)k 


-E 


(PaOp) 


k= 1 


(48) 

(49) 


where  the  last  equality  follows  from  the  Taylor  series  expansion  of  log(l  —  x )  for  x  <  1.  From 
equation  49,  we  have 


lim  nD(B0\\B1) 

n— >oo 


a 


Act  log - b  Api  —  Act 

Pi 

D(Po\\Pi). 


(50) 

(51) 


5.  Confidence  Interval  of  ROC 


In  this  section,  we  provide  a  confidence  interval  of  the  receiver  operating  characteristic  (ROC). 
The  confidence  interval  captures  the  variability  of  the  ROC  across  realizations  of  the  sensor 
network.  The  confidence  interval  of  the  ROC  is  derived  using  Holder’s  inequality.  It  is  a  bound 
and  not  necessarily  tight,  but  applies  in  very  general  cases  where  the  uncertainty  in  the  false  alarm 
probability  (i.e.,  horizontal  variation  in  the  ROC)  is  insignificant  compared  that  of  the  detection 
probability. 

Here,  the  sensor  observations  are  independent  of  the  target  under  the  null  hypothesis,  thus  there  is 
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no  variation  in  the  global  false  alarm  probability  due  to  a  random  target.  This  is  also  true  relative 
to  the  randomness  of  sensor  node  locations  under  the  null  hypothesis.  The  remaining  variability 
in  the  global  false  alarm  probability  is  due  to  the  randomness  in  the  number  of  active  nodes.  For 
the  confidence  interval,  we  assume  there  is  no  horizontal  variation  in  the  ROC  due  to  the 
randomness  in  the  number  of  active  nodes.  For  a  given  number  of  active  nodes,  the  count 
statistics  at  the  fusion  node  are  binomial  (under  perfect  communications).  Since  the  average 
number  of  false  alarms  from  the  sensor  nodes  is  held  constant  (i.e.,  a  — *  0  as  A  — >  oo  s.t.  A  a  is 
constant),  these  binomial  distributions  are  well  approximated  by  Poisson  distributions  as  the 
density  of  the  network  increases.  Thus,  the  confidence  interval  given  here  is  a  measure  of  only 
the  vertical  variance  of  the  ROC. 


For  random  variable  W  G  W  and  mapping  g  :  W  — >  [a,  b],  with  —  oo  <  a  <  b  <  oo,  the  variance 
of  g(W)  can  be  upper  bounded  by 


var (g(W))  <  E(g(W)) 


lim  (  [  (g(w))qdFw) 

\Jw  J 


E  (g(W)) 


(52) 


Let  P(t\  represent  the  (random)  probability  of  detection  for  decision  rule  S(t).  The  randomness  in 
P(i)  follows  from  variability  in  the  number  of  active  nodes,  the  locations  of  the  sensor  nodes  and 
the  signal  source,  and  the  signal  intensity.  Consider  P(l}  as  a  mapping  of  a  random  vector  to  the 
interval  [0,1].  The  variance  of  the  probability  of  detection  is  upper  bounded  by 


var(Pw)  <  /%  (1  -  /%)  . 


(53) 


As  the  probability  of  missed  detection  tends  toward  the  extremes,  the  confidence  interval  on  the 
detection  probability  shrinks.  Figure  3  plots  an  hypothetical  ROC  curve  with  corresponding 
confidence  interval  derived  from  equation  53.  For  the  example  in  figure  3  and  assuming  a  beta 
probability  distribution  for  P(l}  with  mean  (1  —  B{i))  and  variance  B(i){  1  —  /%)),  more  than  99.8% 
of  the  outcomes  have  a  detection  probability  greater  than  0.5  with  a  false  alarm  probability  of 
2  x  10-5  (the  average  detection  probability  being  approximately  0.88). 
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Figure  3.  Notional  ROC  with  lcr  confidence  bounds. 


6.  Numerical  Studies 


In  this  section,  we  demonstrate  system  performance  through  numerical  examples.  The  examples 
exhibit  performance  in  terms  of  ROC  curves,  energy  consumption,  and  probability  of  detection 
versus  sensor  network  density  and  number  of  communicating  slots. 

The  simulations  parameters  for  the  numerical  example  are  as  follows.  The  signal  attenuation  in 
the  propagation  medium  is  given  by 


9(d)  =  ( d 2  +  l)"1  (54) 

where  d  is  the  distance  between  the  sensor  and  signal  source.  This  signal  propagation  loss  model 
is  very  similar  to  those  used  in  references  1,  2,  4,  5.  The  sensor  observations  under  ho  are 
chi-square  distributed  with  10  degrees  of  freedom.  Under  hi,  the  sensor  observations  are 
(conditionally)  noncentral  chi-square  distributed  with  10  degrees  of  freedom  and  noncentrality 
parameter  lO(Asg  (D))2,  given  the  signal-sensor  distance  D  and  signal  amplitude  As.  The  sensor 
locations  are  uniformly  distributed  on  a  disk  of  radius  r  =  20  units.  The  range  and  amplitude  of 
the  signal  source  are  modeled  as  discrete  random  variables.  The  distributions  of  the  signal 
source’s  range  Rs  from  the  center  of  the  disk  and  amplitude  As  are  discrete  uniforms  on  the 
intervals  [0,  20]  and  [2,  200],  respectively.  The  average  number  of  false  alarms  is  Xa  =  10_1. 
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First,  we  compare  the  analytic  ROC  curves  against  sample  estimates  from  105  Monte  Carlo  trials 
for  both  cases:  perfect  communications  and  the  delay-constrained  MAC  protocol.  The  ROC 
curves  represent  the  global  probability  of  detection  versus  global  probability  of  false  alarm.  The 
purpose  is  to  validate  the  analytic  expressions  since  their  evaluation  involves  numerical 
integration  (due  to  the  expectation  with  respect  to  random  sensor  locations).  Figures  4-6  are  plots 
of  ROC  curves,  each  with  four  different  average  number  of  active  nodes:  A  =  10  (blue),  20 
( meen),  40  (red),  80  (  )  nodes.  Figure  4  contains  ROC  curves  for  the  perfect  communications 

case.  The  solid  lines  represent  the  analytic  curves  and  the  circles  are  sample  estimates  of  the 
ROC  curves.  Figure  5  contains  ROC  curves  for  the  MAC  protocol  with  M  =  3  slots.  The 
dashed  lines  in  figure  5  represent  the  analytic  curves  and  the  circles  are  sample  estimates  of  the 
ROC  curves.  As  seen  in  figures  4  and  5,  the  sample  estimates  align  well  with  the  analytic  forms. 
Each  also  demonstrate  how  the  ROC  curves  improve  with  average  number  of  nodes.  The 
remaining  simulation  examples  are  evaluated  from  the  analytic  forms  only. 

Figure  6  overlays  the  analytic  ROC  curves  of  perfect  communications  (solid  lines)  and  the  MAC 
protocol  (dashed  lines)  with  M  —  3  slots.  The  ROC  curves  in  figure  6  are  plotted  on  a  smaller 
domain  of  probability  of  false  alarm  (i.e.,  a  €  (0, 10  3])  to  show  the  slight  reduction  in 
performance  from  assuming  perfect  communications  to  the  MAC  protocol  case,  despite  having 
only  M  =  3  slots  for  the  protocol.  Figure  9  demonstrates  the  detection  probability  versus  a  wider 
range  of  slots  M.  Recall  that  the  fusion  rules  for  each  case  are  different  and  the  fusion  rule  for 
the  MAC  protocol  attempts  retain  as  much  information  as  possible  from  packet  collisions. 

Fusion  rules  in  the  literature  typically  disregard  slots  with  collisions. 

Figure  7  represents  the  normalized  per-node  energy  as  a  function  of  the  average  size  of  the  sensor 
network.  For  this,  we  exclude  the  energy  consumed  by  the  fusion  node  since  that  amount 
remains  constant.  In  particular,  the  normalized  per-node  energy  is  given  by 

=  1  +  Poa  +  (1  —  Po)/3.  (55) 

We  do  not  assume  a  prior  probability  on  the  presence  of  a  signal  source.  Thus,  plotted  in  figure  7 
are  the  extreme  cases  with  p0  =  0  (solid  blue)  and  p0  =  1  (dashed  green).  For  either  case,  or  any 
Po  £  (0,1),  this  shows  that  the  per-node  energy  expenditure  decreases  as  the  average  number  of 
sensor  nodes  increases.  In  particular,  this  is  due  to  desensitizing  the  sensor  nodes  and  therefore 
fewer  nodes  report  detections  to  the  fusion  node. 

Figure  8  demonstrates  the  improvement  in  the  probability  of  detection  as  the  average  number  of 
active  nodes  increases.  Plotted  are  the  probability  of  detection  for  perfect  communications  (solid 
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line)  and  the  MAC  protocol  (dashed  lines)  for  M  —  3  (  ),  5  (o),  and  7  (  )  slots,  and  with  a 
probability  of  false  alarm  a  =  10-8.  The  probability  of  detection  curve  appears  to  flatten 
between  1000  and  3000  nodes.  Although  not  shown,  this  is  due  to  the  choice  of  (discrete) 
distributions  on  the  signal  source.  Roughly,  about  90%  of  the  signal  sources  have  a  probability  of 
detection  that  increases  from  nearly  0  to  1  between  10  and  1000  nodes,  while  the  other  10%  do 
the  same  but  for  more  than  3000  nodes.  However,  according  to  section  4.2.1,  the  increase  in  the 
detection  probability  is  still  strict  since  the  local  ROC  curves  are  differentiable  for  this  simulation 
example.  It  would  not  be  a  surprise  that  the  support  of  and  probability  distribution  of  the  signal 
source  would  have  a  significant  impact  on  the  performance  of  such  a  system. 

Finally,  figure  9  shows  the  probability  of  detection  as  a  function  of  the  number  of 
communications  slots  M  for  the  MAC  protocol  (squares)  for  A  =  10  (blue),  20  (green),  40  (red), 
80  (  )  nodes.  The  probability  of  false  alarm  is  a  =  10-8.  The  solid  lines  represent  the 

probability  of  detection  assuming  perfect  communications,  which  does  not  depend  on  the  number 
of  slots.  Figure  9  demonstrates  that,  even  in  the  presence  of  collisions,  performance  approaches 
that  of  the  ideal  communications  case  for  increasing  number  of  slots.  However,  increasing  the 
number  of  slots  also  increases  the  delay  at  the  fusion  node  to  make  a  global  decision. 
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Figure  4.  ROC  curves  for  A  =  10  (blue),  20  (  reen),  40  (red), 
80  (  )  for  perfect  communications.  Solid  lines 

represent  the  analytic  ROC  and  circles  represent 
sample  estimates. 


Figure  5.  ROC  curves  for  A  =  10  (blue),  20  (green),  40  (red), 
80  (  )  for  the  MAC  protocol  with  M  =  3  slots. 

Dashed  lines  represent  the  analytic  ROC  and  circles 
represent  sample  estimates. 
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Figure  6.  ROC  curves  for  A  =  10  (blue),  20  (green),  40  (red), 
80  (  ).  Solid  and  dashed  lines  represent  the 

analytic  ROC  for  perfect  communications  and  the 
MAC  protocol  with  M  =  3  slots,  respectively. 


Figure  7.  Normalized  average  per-node  energy  (Es/ (A E)) 
versus  number  of  active  nodes  A  assuming  7To  =  0 
(solid  blue)  and  ttq  =  1  (dashed  green). 
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Figure  8.  Global  probability  of  detection  versus  number  of 

active  nodes  A  for  Pfa  =  10-8.  The  curves  represent 
Pd  for  perfect  communications  (blue)  and  the  MAC 
protocol  with  M  =  3  (green),  M  =  5  (red),  and 
M  =  7  (  )  slots. 
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Figure  9.  Global  probability  of  detection  versus  number  of 

communications  slots  M  for  Pfa  =  10-8  and  A  =  10 
(blue),  20  (green),  40  (red),  80  (  ). 
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7.  Conclusion 


We  derived  the  fusion  rules  for  the  case  of  sensor  nodes  communicating  local  binary  decisions  to 
a  common  fusion  node  for  (1)  assuming  perfect  communications  and  (2)  through  a 
delay-constrained  random-access  MAC  on  a  single-hop  wireless  network.  The  fusion  node 
combines  received  information  to  make  a  global  decision  on  the  presence  or  absence  of  a  signal 
source.  Under  the  MAC  protocol,  there  is  a  chance  that  some  detection  messages  are  not  received 
properly  due  to  collisions  with  other  messages  since  the  communications  medium  is  shared.  In 
this  work,  sensor  nodes  attempt  to  transmit  messages  representing  local  detections  only.  Thus, 
collisions  contain  useful  information  about  local  detections.  The  fusion  rules  simply  rely  on 
count  statistics:  (1)  the  number  of  detections  for  perfect  communications  and  (2)  the  number  of 
successful  transmissions  and  the  number  of  collisions  for  the  MAC  protocol.  The  perfect 
communications  case  is  seen  as  an  asymptote  of  the  MAC  protocol  case  (i.e.,  as  M  — y  oo). 

The  local  decision  rules  are  constrained  to  be  identical.  Additionally,  the  system  is  constrained  to 
maintain  a  constant  number  of  false  alarms  across  the  sensor  network  in  any  given  activity  period. 
This  implies  that  the  sensor  nodes  are  desensitized  as  the  average  number  of  active  nodes 
increases.  We  showed  that,  for  local  decision  rules  that  have  concave  ROC  curves,  the  global 
detection  performance  improves  with  average  number  of  active  nodes,  despite  suppressing  local 
decisions. 

It  was  shown  that  the  performance  under  the  MAC  protocol  is  somewhat  degraded  relative  to  the 
case  of  “ideal"  communications.  However,  simulations  show  that  the  gap  in  performance  can  be 
negligible,  despite  collisions  (without  retransmissions)  in  communicating  detection  messages. 
Additionally,  the  gap  in  performance  can  be  reduced  by  increasing  the  number  of  slots  M,  but  at 
the  expense  of  added  delay. 

In  addition,  we  provided  a  bound  on  variance  of  the  ROC  curves.  We  assumed  a  random  sensor 
network  and  signal  source.  The  ROC  curve  for  a  given  sensor  network  and  signal  source  will 
vary  from  the  average  case.  While  the  bound  can  be  loose,  it  does  not  depend  on  the  actual 
distributions  and  is  simple  to  calculate  for  a  given  average-case  ROC  curve. 

This  framework  could  be  extended  to  an  asynchronous  (or  pure)  ALOHA  to  reduce  energy 
consumption  in  the  sensor  network.  S-ALOHA  requires  time  synchronization  between  all  the 
sensor  nodes  and  the  fusion  node.  Additional  procedures  and  communications  are  needed  beyond 
those  for  signal  detection  to  maintain  synchronization  within  a  desired  tolerance.  Removing  the 
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need  for  synchronization  also  eliminates  an  energy  budget  associated  with  synchronization. 
However,  it  is  well  known  that  S-ALOHA  reduces  the  probability  of  a  collision  over  pure 
ALOHA  due  to  halving  of  the  vulnerability  period.  Thus,  it  is  reasonable  to  conclude  that  the 
detection  performance  would  be  degraded  compared  to  the  MAC  protocol  discussed  here. 
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